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We calculate the spin susceptibility in the s± and s++ superconducting states of the iron pnictides 
using the effective five orbital model and considering the quasiparticle damping. For the experimen- 
tally evaluated magnitude of the quasiparticle damping and the superconducting gap, the results 
at the wave vector ~ (tt, 0) show that the s± state is more consistent with the neutron scattering 
experiments, while for larger quasiparticle damping and the superconducting gap, the s++ state can 
be more consistent. To distinguish between two cases that reproduce the experiments at the wave 
vector ~ (tt, 0), we propose to investigate experimentally the wave vector ~ (tt, tt). 



PACS numbers: 74.20. Rp, 78.70.Nx, 74.70.Xa 

The discovery of the iron-based superconductors has 
received considerable attention. The high transition tem- 
perature Tjji itself has given a dramatic impact, but the 
possibility of peculiar unconventional pairing state has 
also been an issue of great interest. In fact, it has been 
proposed theoretically at the early staged— , that a nat- 
ural pairing state is the so called s±-wave pairing, where 
the superconducting gap is fully open, but changes its 
sign across the wave vector that bridges the disconnected 
Fermi surfaces. The sign change occurs because the spin 
fluctuations that develop at the Fermi surface nesting 
vector gives repulsive pairing interaction. 

There are many experimental results which suggest 
that the order parameter is fully gapped in a number of 
iron-based materials, such as the penetration depth 8 - - — 
and the angular-resolved photo-emission spectroscopy 
(ARPES).— ~— There are also experimental suggestions 
that the iron-based materials are unconventional super- 
conductors. The nuclear magnetic relaxation rate lacks 
coherence peak below Tj£~—, which suggests the pres- 
ence of sign change in the superconducting gap. Also in 
ref integer and half-integer flux-quantum transitions 
in composite niobium-ion pnictide loops have been ob- 
served, which suggests the presence of the sign change in 
the superconducting gap. In refi±2, STM/STS measure- 
ments have been performed to detect the quasiparticle in- 
terference indicating the realization of the s± gap. These 
experiments seem to be consistent with the s± scenario. 
However, recently there has been some debate concerning 
the sensitivity of T c against impurities. Although some 
experiments show suppression of T c by impurities^, it 
has been pointed out, e.g., in refs. 21 ' 22 that the suppres- 
sion of T c by impurities is too weak for a pairing state 
with a sign change in the gap. A calculation based on 
a five band model by Onari et al. has supported this 
theoretically 2 ^, although the strength of the impurity po- 
tential adopted there is large compared to those calcu- 
lated from first principles 24 . As a possible pairing state 
that is robust against impurities, the so-called s ++ state, 



where the gap does not change its sign between the Fermi 
surfaces, has been proposed^—. 

Now, it has been proposed at the early stag o 5 ' 28 ' 29 that 
one of the promising ways to determine whether the su- 
perconducting gap indeed changes its sign between the 
disconnected Fermi surfaces is the observation of neutron 
scattering resonance at the nesting vector of the elec- 
tron and hole Fermi surfaces. In fact, neutron scattering 
experiments have indeed observed a peak like structure 
in the superconducting stated— . This has been taken 
as a strong evidence for the sign change in the super- 
conducting gap. However, Onari et al. later took into 
account the quasiparticle damping effect in the calcu- 
lation of the dynamical spin susceptibility, and showed 
that a peak like enhancement over the normal state val- 
ues can be seen even in the s++ state, which is due to 
the suppression of the normal state susceptibility orig- 
inating from the damping 3 -^ In ref.—, the strength of 
the quasiparticle damping was estimated from the ex- 
perimental results to be ~ 10(meV) or less, but larger 
values (> 50(meV)) was adopted in the actual calcula- 
tion, fixing the gap/damping ratio to be around unity, 
due to the restriction in the numerical calculation. 

In the present study, we revisit the problem of the res- 
onance peak in the neutron scattering experiment. We 
adopt the same formalism as in refj 3 ^, but with smaller 
and realistic values (~ 10(meV)) of the quasiparticle 
damping and the superconducting gap by taking fe-point 
meshes up to 16384x 16384. For such small values of 
the quasiparticle damping and the superconducting gap, 
the resonance peak enhancement over the normal state 
susceptibility at the wave vector ~ (71", 0) is found to be 
comparable to those observed experimentally, while the 
enhancement over the normal state susceptibility for the 
S4.4. state remains to be small. As we increase the quasi- 
particle damping and the superconducting gap up to val- 
ues of an order of magnitude larger, the enhancement 
over the normal state susceptibility comparable to the 
experimental results arises in the s++ state, in agreement 
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FIG. 1: (Color online) Energy dependence of x"(Qi E) at Q = 
(it, 7t/16) with the quasiparticle damping of 70 = 10(meV). s± 
(circle), s++ (triangle), and normal (square) states. We take 
(a)A = 5(meV) and U = 1.375(eV), and (b)A = 10(meV) 
and U = 1.3(eV). 



with the results in refj^ 5 -. For such a large quasiparticle 
damping and superconducting gap, an enhancement also 
occurs at the wave vector (tt, it). This is in contrast to the 
case with the small damping, where such an enhancement 
does not occur at (71", tt) (regardless of the pairing state). 
Thus, we propose that by looking also at the wave vector 
(•7r, 7r) in the neutron scattering experiments, one can dis- 
tinguish the origin of the "resonance like" enhancement 
at (tt,0). 

We apply multi-orbital RPA (Random phase approxi- 
mation) to the five orbital model of LaFeAsO obtained in 
the unfolded Brillouin zonei, where the x- and y- axes are 
taken in the Fe-Fe bond direction. The orbital-dependent 
spin susceptibility x(Q,E) is given as 



C (Q, E)=T f(i - Xo(Q, Efi^xaiQ, E)] " , (1) 

L — ' L J cc 



where the interaction coefficients U c ^ hcd — U, U', J, and 
J' for a = b = c = d 1 a = c=/=b = d, a = b =/= c = d, 
and a = d 7^ b = c, respectively. Throughout this paper, 
we put J = J' = 0.15(eV), U 1 = U - 2J and fix the 
chemical potential at [i = 10.93(eV) , which corresponds 
to an optimally doped regime of n ~ 6.1 (10% doping). 



Xo(Q,E) is the bare spin susceptibility expressed as 
[Xo(Q, E)C = -J2Yl [ M 2'cf{k, k + Q)xZik, k + Q,E) 

k vv ! 

+M%g(k, k + Q) X v i (k, k + Q, E) . 

(2) 

Here, XoG(F)(k,k + Q,E) denotes the normal (anoma- 
lous) part of the band-dependent BCS spin susceptibility 
written as 



„f|2|„,i/' 



XZ (k,k + Q,E) = 



X7 F (k,k + Q,E) 



(3) 
(4) 



at zero-temperature (E > 0) with E v k = ^ef + |A£| 2 , 



,i/|2 



(1 + el/El)/2, |^| 2 = (1 - ?JEl)/2, ulvl* 
A^/(2_E^), and e£ is the v-th. band energy measured rela- 
tive to the Fermi energy. M^' c f{k,k+Q) {M™'J(k,k + 
Q)) is given by 



M. 



X' d G (fc, k + Q) = U* av {k)U bvl {k + Q)U; v ,{k + Q)U du {k), 

(5) 



M y v 'F 



abcd (fe, k + Q) = u:„{k)U bv ,{k + Q)U c „(k)U* d „,(k + Q), 

(6) 



with the unitary matrix U (k) which diagonalizes the 
Hamiltonian in the orbital basis. Here, we introduce 
the band-index v whose energy t v satisfies the relation 
e v > e v (y > v'). There are the hole Fermi surfaces on 
the 2nd and 3rd bands around (k x ,k y ) = (0,0) and the 
electron Fermi surfaces on the 4th band around (tt, 0) and 
(tt, 0). For the 's ++ -wave', we take A 2 = A 3 = A 4 = A 
and for 's±-wave' A 2 = A 3 = —A 4 = A . As done 
in ref.— , we introduce a Gaussian cutoff for the gap 
A£ = A^exp{-[e£/A£;] 2 }, and take AE = 4A . 

The quasiparticle damping is taken into account in the 
same form as in ref. 35 : T W M q = max {j(E l k ),j(E l k+q )} 
with 7(e) = a(e)7o where 7(e) = r\ for |e| < 3Ao, 7(e) = 
70 for |e| > 4A , and linear interpolation for 3A < |e| < 
4Ao^ In the present study, the smearing factor is taken 
to be 77 = 0.5(meV) in order to consider small values 
(~ 10(meV)) of 70. To cope with the realistic magnitude 
of the superconducting gap, we take 8192 x 8192 fc-point 
meshes throughout the paper. We have confirmed that 
the results do not change if we take 16384 x 16384 k- 
point meshes. 

We first calculate x"{Q-,E) = lmx(Q,E) at the nest- 
ing vector of the electron and hole Fermi surfaces Q = 
(71,71-/16) for s± and s ++ states. We adopt the values 
A = 5(meV) in (a) and A = 10(meV) in (b) for the 
superconducting gap. We take 70 = 10(meV) because 
70 is estimated to be of the same order as Arj 3 - 5 . The 
values of U are determined so that the normal state sus- 
ceptibility is broadly maximized around several times of 
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Ao as observed experimentally. As shown in Fig. [T] the 
resonance peak develops above the normal state value for 
the s±-wave state, which is qualitatively consistent with 
the previous studies 5 - ' 28 ! 29 . In contrast to the case with 
larger quasiparticle damping and superconducting gap 35 , 
the spectral weight at the peak is only about three times 
larger compared to the normal state value, which can 
be considered as consistent with the experimental results 
^11 J±. We note here that not only the small Ao and 70, 
but also the Gaussian cutoff for the superconducting gap 
is also effective in reducing the spectral weight in the su- 
perconducting stated. On the other hand, there is only 
a small enhancement in the case of s++ again in contrast 
to the case with larger quasiparticle damping and super- 
conducting gap- 5 . Thus, within the present formalism, 
as far as we adopt the 70 value estimated from experi- 
mental results^ 5 - and the realistic magnitude of the su- 
perconducting gap, the s± state is more consistent with 
the experiments. 

To see how the peak like structure develops in the s++ 
state for stronger quasiparticle damping, we now calcu- 
late X "(Q = 0", 71-/16), -E) with 7o = 2A = 50(meV) as 
shown in Fig. [Ha). We see that the peak structure in- 
deed develops for such a large 70 and Ao, consistent with 
ref.—. The origin of the enhancement above the normal 
state values is mainly due to the the strong suppression 
of x" m the normal state by the quasiparticle damping. 

For such a large 70 and Ao, the enhancement ratio 
of the resonance peak in the s± state against the nor- 
mal state value becomes somewhat larger as shown in 
Fig. [Hh). In fact, for even larger values of Ao and 70, 
the resonance peak in the s± state becomes much larger 
compared to the experiments, as was shown in ref. 35 . 
Thus, as we increase the values of 70 and Ao by an order 
of magnitude larger than the experimentally evaluated 
values, there is a tendency that the s++ state becomes 
more consistent with the experiments than the s± state. 

So far, our results for the wave vector Q = (it, 7r/16) 
show that the s± state is consistent with the experiments 
when small values of 70 and Ao are adopted, while s++ 
can explain the experiments for large 70 and Ao. The 
former parameter values are realistic as far as the ex- 
perimental evaluations are concerned, but there might 
remain a possibility that the effective values of 70 and/or 
Ao that should be adopted in the present formalism 
turn out to be larger—. In that sense, the possibility 
of s ++ state may not be ruled out. Here we propose 
a way to distinguish whether the peak like structure 
at ~ (tt, 0) originates from the resonance effect in the 
s± state or the strong quasiparticle damping effect in 
the s ++ state. We calculate x"(Q>E) at Q = (ir,n) 

Q = (tt, tt) is the translation vector for which the 
two electron Fermi surfaces around (it, 0) is nearly su- 
perposed to that around (0,7r). (It should be noted 
however that this wave vector is not the nesting vector 
that induces spin fluctuations since this translation su- 
perposes an electron Fermi surface on another electron 
Fermi surface.) We consider the two cases where the ex- 
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FIG. 2: (Color online) Energy dependence of x"(Q,E) at 
Q — (7r, 7r/16) in the case of (a) s++ and (b) s± with Ao = 
25(meV), the quasiparticle damping 70 = 50(meV), and U = 
1.3(eV). 



perimental results can be explained at Q = (tt, 7t/16) 
: (a) s±-wave with 7 = 2A = 10(meV) and (b) s++- 
wave with 7 = 2Ao = 50(meV). Here we consider only 
the case of s± (s++) for small (large) 7 and A because 
this is the combination that is consistent with the experi- 
ment at ~ (tt, 0), but we note here that actually, for each 
fixed parameter set, the calculation results at (tt,tt) are 
barely affected by whether the gap is s± or s++ because 
in both cases the superconducting gap does not change 
its sign between the Fermi surfaces separated by (tt,tt). 
As shown in Fig. [21a) , no enhancement above the nor- 
mal state values is found in the former case because the 
quasiparticle damping is small. On the other hand, in the 
latter case, we find a peak structure developing above the 
normal state values even at (tt, tt) because the effect of 
the quasiparticle damping is strong. Therefore, by inves- 
tigating the spin fluctuations at ~ (tt, it) using neutron 
scattering experiments, one can distinguish the origin of 
the peak like structure at ~ (tt, 0) ; if an enhancement 
above the normal state values is found in the supercon- 
ducting state at (tt, tt) as well as at (tt, 0), its origin may 
well be due to the quasiparticle damping, while if an en- 
hancement is observed at (tt,0) but not at (tt, tt), the 
quasiparticle damping is not so strong and the origin is 
the resonance in the s± state. 

In conclusion, we have calculated the spin susceptibil- 
ity in the s± and s++ superconducting states of the iron 
pnictides by applying multi-orbital RPA to the effective 
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FIG. 3: (Color online) Energy dependence of x"{Q,E) at 
Q = (tt,tv) for parameter values where (a)s± and (b)s++ 
states are consistent with the experiments at the wave vector 
~ (tt,0). (a)A = 5(meV), U = 1.375(eV) and 7o = 10(meV), 
and (b)A = 25(meV), U = 1.3(eV) and 7o = 50(meV). 



five-band model and considering the quasiparticle damp- 
ing. We have found that as far as we adopt the values of 
the quasiparticle damping and the superconducting gap 
evaluated from experiments, the enhancement above the 
normal state values at ~ (n, 0) observed in the neutron 
scattering experiments is more consistent with the reso- 
nance peak in the s± state. On the other hand, for large 
magnitude of the quasiparticle damping and the super- 
conducting gap, the enhancement in the spin suscepti- 
bility of the state over the normal state values can 
explain the experimental results^, while the resonance 
peak in the s± state becomes large compared to exper- 
imental observations. In case there is a possibility that 
the quasiparticle damping and the superconducting gap 
that should be adopted in the present formalism is effec- 
tively larger than the experimentally evaluated values, 
we propose to investigate in the neutron scattering ex- 
periments the wave vector around (tt, it) in the unfolded 
Brillouin zone to see if there is an enhancement like those 
found at (it, 0) in the superconducting state. 
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The present calculation has been done for a purely two 
dimensional model of LaFeAsO. On the other hand, the 
122 materials exhibit strong three dimensionality, and a 
possibility of superconducting gap nodes or dips on the 
strongly warped portion of the hole Fermi surface has been 
pointed out theoretically ( S. Graser et al, Phys. Rev. B 
81, 214503 (2010) and K. Suzuki et al, J. Phys. Soc. Jpn. 
80, 103710 (2011)) or experimentally (e.g., F. Hardy, et 
al, EuroPhys. Lett. 91 (2010) 47008). In ref.^i, we have 
performed a calculation for a three dimensional model of 
122 material (without considering quasi particle damping), 
where taking into account the above mentioned effect tends 
to reduce the strength of the resonance peak in the s± 
state, which gives further agreement with experiments. 
For instance, the energy scale of the superconducting tran- 
sition temperature and the superconducting gap in the ran- 
dom phase approximation tends to be larger compared to 
actual experimental values. 



